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Abstract 

Background: Colorectal adenoma develops into cancer with the accumulation of genetic and epigenetic changes. We 
studied the underlying molecular and clinicopathological features to better understand the heterogeneity of colorectal 
neoplasms (CRNs). 

Methods: Vile evaluated both genetic (mutations of KR AS, BRAF, TP53, and PIK3CA, and microsatellite instability [MSI]) and 
epigenetic (methylation status of nine genes or sequences, including the CpG island methylator phenotype [CIMP] markers) 
alterations in 1 58 CRNs including 56 polypoid neoplasms (PNs), 25 granular type laterally spreading tumors (LST-Gs), 48 non- 
granular type LSTs (LST-NGs), 19 depressed neoplasms (DNs) and 10 small flat-elevated neoplasms (S-FNs) on the basis of 
macroscopic appearance. 

Results: S-FNs showed few molecular changes except SFRP) methylation. Significant differences in the frequency of KRAS 
mutations were observed among subtypes (68% for LST-Gs, 36% for PNs, 16% for DNs and 6% for LST-NGs) (P<0.001). By 
contrast, the frequency of TP53 mutation was higher in DNs than PNs or LST-Gs (32% vs. 5% or 0%, respectively) (P<0.007). 
We also observed significant differences in the frequency of CIMP between LST-Gs and LST-NGs or PNs (32% vs. 6% or 5%, 
respectively) (P<0.005). Moreover, the methylation level of LINE-1 was significantly lower in DNs or LST-Gs than in PNs 
(58.3% or 60.5% vs. 63.2%, P<0.05). PIK3CA mutations were detected only in LSTs. Finally, multivariate analyses showed that 
macroscopic morphologies were significantly associated with an increased risk of molecular changes (PN or LST-G for KRAS 
mutation, odds ratio [OR] 9.1 1; LST-NG or DN for TP53 mutation, OR 5.30; LST-G for PIK3CA mutation, OR 26.53; LST-G or DN 
for LINE-1 hypomethylation, OR 3.41). 

Conclusion:\Ne demonstrated that CRNs could be classified into five macroscopic subtypes according to clinicopathological 
and molecular differences, suggesting that different mechanisms are involved in the pathogenesis of colorectal 
tumorigenesis. 
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Introduction 

Colorectal cancer (CRC) can develop via various molecular 
pathways. Most CRCs develop over a long period of time through 
a multistep process called the adenoma-carcinoma sequence [1]. 
Approximately two-thirds of sporadic CRCs arise from conven- 
tional adenomas and usually show a protruding (polypoid) 
macroscopic appearance. The process of colorectal carcinogenesis 
often begins with the inactivation of the ^PC/fi-catenin signaling 
pathway (the Vogelstein model), followed by KRAS and TP53 
mutations [2]. However, serrated adenomas (SAs), particularly 



sessile serrated adenoma/polyias (SSA/Ps), have been described as 
the immediate precursors for CRCs that develop via an alternative 
pathway consisting of the CpG island methylator phenotype 
(CIMP) and BRAF mutations [3,4] . In addition, CIMP cancers 
may develop either via a mutator (microsatellite instability; MSI) 
pathway, or via a pathway that leads to microsatellite stability 
(MSS) [5]. However, additional pathways that are not fully 
vmderstood may also contribute to colorectal carcinogenesis. 

Recent studies [6-8] revealed that conventional (non-serrated) 
adenomas could be morphologically classified into polypoid 
neoplasms (PNs) and nonpolypoid neoplasms (NPNs; also referred 
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to as flat and depressed neoplasms). PNs develop via the traditional 
adenoma-carcinoma sequence, and their tumorigenesis is charac- 
terized by loss of heterozygosity, which leads to the inactivation of 
tumor suppressor genes such as APC and TP 5 3 [2] . The mutation 
of KRAS, which activates the mitogen-activated protein kinase 
(MAPK) cascade and promotes malignant transformation, is a key 
event in the adenoma-carcinoma sequence [9] . In contrast, NPNs 
have a low frequency of KRAS mutation, and the chromosomal 
changes that occur during the development of NPNs are markedly 
different from those that occur during the progression of PNs [10- 
13]. We prc'\'i()usly showed that NPNs have a higher frequency of 
MSI (MSI-H), an increased abnormal accumulation of phosphor- 
ylated MAPK protein, and a lower frequency otKRAS mutations 
than PNs [14]. In addition. An et al. [15] reported specific patterns 
of aberrant DNA methylation in CIMP-negative CRCs, particu- 
larly a decrease in global DNA methylation, and an increase in the 
age-related methylation of multiple genes such as MGMT, 
RASSFl and SFRPl. 

Clinically, depressed neoplasms (DNs) are characterized by 
increased risk of malignancy as compared with PNs, even if they 
are small [16]. In addition, flat cl(^\-atcd neoplasms can be 
classified into small-flat adenoma and latc'rally spreading tumors 
(LSTs), which were initially reported by Kudo el al. [7], and are 
characterized by lateral extension along the luminal wall with a 
low vertical axis. These tumors are sub-categorized into granular 
type LST (LST-G) and non-granular t\ipe LST (LST-NG), based 
on different molecular features [17,18]. NPNs are not easily 
detected during a colonoscopy and it is challenging for the 
colonoscopists to distinguish them from the normal mucosa [7]. 
Inadequate recognition of NPNs could result in the development 
of interval cancers. Molecular analysis of precursor lesions and 
early stage CRC should therefore be performed to gain a better 
understanding of the different pathways of CRC de\ (^lopment. In 
addition, such analyses may be used for the implementation of the 
appropriate screening and therapeutic intervention programs for 
CRCs. 

We hypothesized that the epigenetic and genetic features of 
CRC may be shared with macroscopic subtypes of conventional 
neoplasms such as adenomas and Tl carcinoma. In the present 
study, we tested this hypothesis by determining the genetic and 
epigenetic profiles of colorectal neoplasms (CRNs) by using 
molecular markers associated with the mechanism of colorectal 
carcinogenesis. 

Materials and Methods 

Patients and samples 

A total of 158 CRNs from 153 patients who underwent 
endoscopic (n=125) or surgical resection (n = 33) at Showa 
University Hospital were examined. The samples were selected 
solely based on tissue availability. Patients with inflammatory 
bowel disease or with a familial predisposition to cancers such as 
familial adenomatous polyposis or hereditary nonpol^'posis colo- 
rectal cancer were excluded. The ethics committee of the Showa 
University School of Medicine approved the procedures for tissue 
collection and analysis, and written informed consent was obtained 
from each patient. 

Endoscopic evaluation and macroscopic classification 

AU patients were prepared for the procedure by administration 
of 1.8 L oral electrolyte lavage solution. Colonoscopists with 
extensive experience performed all examinations by using high- 
resolution video colonoscopes (CF-240ZI, CF-260AI, or CF- 
260HZI; Olympus Optical Co., Tokyo, Japan). CRNs were then 



prospectively classified as PN (Figure SI a) or NPN (Figure Slb-e) 
based on the Paris endoscopic classification [19]. Briefly, NPNs 
were defined as neoplasms showing slight mucosal elevation with a 
flat or slightly rounded surface and a height of less than half the 
diameter of the lesion. Histological examination indicated that 
NPNs t)'j)i( ally displayed dysplastic mucosal thickness less than 
twice that of the adjacent nondysplastic mucosa. NPNs comprised 
flat-elevated (Figure Sib, d, e) and DNs (Figure Sic), depending 
on the presence of a depressed component. The flat-elevated 
lesions were subclassified into small flat-elevated neoplasms (S- 
FNs) (Figure Sib), LST-G (Figure Sid) and LST-NG (Figure Sle). 
Briefly, LSTs are defined by a large lateral diameter (>10 mm), a 
low vertical axis, and lateral extension along the luminal wall [20] . 
LST-Gs are composed of superficial spreading aggregates of 
nodules that form flat, broad-based lesions with a granulonodular 
and uneven surface, whereas LST-NGs have a flat smooth surface 
without granulonodular formation. In contrast, PNs presented 
with sessile, pedunculated or semipedunculated morphology 
(Figure SI a). The details of this macroscopic classification based 
on colonoscopic findings are summarized in Table SI in File SI. 

Tissue samples and histological evaluation 

Serial sections (3 \Lm) were obtained from paraffin-embedded 
blocks and prepared for hematoxyhn and eosin (H&E) staining. All 
H&E-stained slides were analyzed by a senior pathologist (T.Y.) 
who was bhnded to the endoscopic findings. Hyperplastic polyps 
and SAs were not included in this analysis. 

To extract genomic DNA, 1 5 formalin-fixed paraffin-embedded 
samples and 143 frozen tissue samples were used. The frozen tissue 
samples were obtained from colonoscopic biopsy specimens and 
stored at -80°C. We distinguished between neoplastic and non- 
neoplastic areas of the specimens based on pit patterns observed 
during chromoendoscopic examination [21,22]. DNA was ex- 
tracted from frozen tissue samples by using standard proteinase 
K/ phenol/ chloroform methods. 

Serial slides were obtained from the archival formalin-frKcd, 
paraffin-embedded tumor tissues. One slide was stained with H&E 
for microdissection. After microdissection, DNA was extracted 
using the QIAamp DNA mini kit (QIAGEN Inc., Valencia, CA). 

Bisulfite polymerase chain reaction and pyrosequencing 
analysis of DNA methylation 

Bisulfite treatment was performed as previously described [23], 
and 2-3 (xL bisulfite-treated DNA was used as the template for 
polymerase chain reaction (PGR). The primers and PCR 
conditions used for the amplification of target genes have been 
described previously. The protocol for pyrosequencing, a quan- 
titative tool for determining methylation density, has been 
described in detail [24]. Pyrosequencing can be used to measure 
the level of methylation at several CpG sites in a given promoter, 
and the methylation status of different sites is usually consistent. 
For each gene, the methylation percentage of all CpGs measured 
was averaged. 

Methylatlon-reiated genes and definition of CIMP 

We studied 9 genes or sequences (MINTl, MINT2, MINT31, 
CDKN2A, MLHl, MGMT, RASSFl, SFRPl, and LINE-1) in 
this analysis. The PCR and sequencing primer sequences used 
were previously reported [24—26]. Sporadic CRCs can be 
classified into 2 groups, CIMP-positive and CIMP-negative 
according to the frequency of methylation of the CpG islands in 
die promoter of 5 genes (MINTl, MINT2, MINTS 1, CDKN2A, 
and MLHl) [27]. As determination of the CIMP status requires a 
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quantitative tool, positive methylation status was defined as a 
methylation density greater than 15% [24]. A tumor was 
considered to be CEVIP-positive if 2 or more GIMP markers were 
methylated, as described previously [27]. All other tumors were 
defined as C IMF-negative. 

KRAS, BRAF, TP53 and PIK3CA gene mutations and MSI 

Samples were analyzed using PCR-based pyro.secjucncing to 
determine thi; pr(\sence of activating mutations in codons 1 2 and 
13 ot: KRAS, codon 600 of BRAF, and in exons 9 and 20 of 
PIK3CA [28-30]. Exons 5 to 8 of TP 5 3 and MSI were assessed 
following previously described protocols [3 1] . 

Statistical analysis 

Pyrosequencing provides a methylation level (%), which was 
analyzed as a continuous variable for the comparison of each gene 
with cUnicopathological variables. The mean, median, ranges, and 
95% confidence intervals (CI) were calculated. 

Differences in the continuous variables (age, tumor size, and 
methylation density) among groups were analyzed by using the 
Kruskal-Wallis test. Post-hoc tests such as the Steel-Dwass method 
were used to compare differences in the continuous variables 
between groups, and P<0.05 was considered to be statistically 
significant. Categorical variables were compared between subtypes 
by using or Fisher's exact test when testing small samples. AU 
tests were two-sided. F values were considered to be significant at a 
Bonferroni-corrected alpha of 0.05. 

Logistic regression analysis using the stepwise method was 
performed to evaluate the relationship between molecular 
alterations of CRNs and gender, age, tumor location (proximal 
vs. distal), tumor size, macroscopic types, and carcinoma 
component (Tl cancer). In this analysis, gender, tumor location, 
macroscopic t^pe, histology, and genetic alterations were consid- 
ered as categorical variables, whereas age and tumor size were 
used as continuous variables. The odds ratio (OR) and 95% CI 
were determined for a variety of factors. P<0.05 was considered to 
be significant. All statistical analyses were performed with SPSS 
version 14.0 (SPSS, Inc., Tokyo, Japan) and JMP version 10 (SAS 
Institute, Inc., Cary, NC). 

Results 

We analyzed the molecular features of 158 CRNs. Table 1 
summarizes the patients' clinicopathological characteristics. No 
significant difierences in the clinicopathological features were 
observed between all CRNs and those with high-grade dysplasia 
(HGD)/submucosal cancer (Tl cancer). 

Clinicopatoioglca! and molecular features of macroscopic 

subtypes of colorectal neoplasms 

Ten out of 158 CRN lesions were S-FNs, and had few 
molecular changes except SFRPl methylation (Table S2 in File 
SI). Therefore, four subtypes were further analyzed. The 
cfinicopatological features of the four phenotypes of CRNs are 
summarized in Table 2. LST-Gs were significandy larger than the 
other subtypes (F<0.05 by Steel-Dwass). Regarding histological 
grade, we observed a significant difference between DNs and the 
other subtypes (F<0.0001). 

The molecular features of four subtypes of CRNs are shown in 
Table 2. Significant dilferences in the frequency of KRAS, TP53, 
PIK3CA mutations and CIMP were observed among the four 
subtypes. We found significant differences in the frequency of 
KRAS mutations among the subtypes (PNs vs. LST-NGs, 
P = 0.0003; LST-Gs vs. LST-NGs, P<0.0001; LST-Gs vs. DNs, 



P = 0.0008; PNs vs. LST-Gs, F = 0.0070, respectively). By contrast, 
the frequency of TPS 3 mutation was higher in DNs than PNs or 
LST-Gs (DNs vs. PNs, P = 0.0066; DNs vs. LST-Gs, P = 0.0038, 
respectively). Significant differences in the frequency of CIMP 
were observed between LST-Gs and LST-NGs or PNs (LST-Gs 
vs. LST-NGs, F = 0.0035; LST-Gs vs. PNs, P = 0.0012, respec- 
tively). Although MGMT and RASSFl methylation density did 
not differ significantly between any subgroups, there were 
difiFerences in the DNA methylation levels of SFRPl between 
two subtypes (P<0.05 by Steel-Dwass except LST-NGs vs. DNs). 
The DNA methylation level of LINE- 1 was significantly lower in 
LST-Gs or DNs than in PNs (P<0.05 by Steel-Dwass). 

PIK3CA mutations in CRNs 

Of 158 CRN samples, 151 were examined for PIK3CA 
mutation status. The frequency of PIK3CA mutation was 
uncommonly low in our studied samples (3%, 5/151; Table 2). 
However, the frequency was higher in LST-G tumors compared 
with the other sub-groups (LST-Gs, 17% [4/23]; LST-NGs, 2% 
[1/45]; PNs, 0% [0/56]; DNs, 0% [0/17]). Of these, die 
difference in the frequency of PIK3CA mutations between LST- 
Gs and PNs was statistically significant (P = 0.0076). Four of the 
five tumors with a PIK3CA mutation were diagnosed as HGD or 
Tl cancer. In addition, PIK3CA mutation was not significandy 
correlated with other molecular changes such as KRAS mutation 
or CIMP. 

Multivariate analysis 

Finally, multivariate analysis was performed to determine 
whether molecular alterations can be predicted by clinicopatho- 
logical variables (Tables 3 and Table S3 in File SI). In the 
multivariate analysis, PN/LST-G morphology and tumor size 
were significant risk factors for KRAS mutation, whereas LST- 
NG/DN morphology was the only significant risk factor for TP33 
mutation. The morphology of LST-G was also the only significant 
risk factor for PIK3CA mutation. Regarding epigenetic alter- 
ations, patient age and tumor size were significant risk factors for 
CIMP. In the multivariate analysis for risk factors for LINE-1 
hypomethylation, LST-G/DN morphology and carcinoma com- 
ponents were significant risk factors. In addition, we validated the 
results of multivariate analysis (Table S4 in File SI). 

Discussion 

In the present study, we identified the distinct features of five 
macroscopic phenotypes of conventional neoplasms based on 
genetic alterations and DNA methylation profiles as follows: PN, 
high frequency of KRAS mutation; LST-G, high frequency of 
CIMP, KRAS and PIK3CA mutation and LINE-1 hypomethyla- 
tion; LST-NG, relatively high frequency of TP 53 mutation; DN, 
high frequency of TP 5 3 mutation and LINE- 1 hypomethylation; 
S-FN, rare methylation changes and genetic alterations. These five 
subtypes are correlated with the diversity of colorectal tumorigen- 
esis. Our data also suggest that the five subtypes of CRNs differ 
clinically. For example, LST-G cases tended to be larger, whereas 
DN lesions more frequently had advanced neoplasms. However, 
S-FNs showed no significant molecular features in our analysis. 
Moreover, whether these macroscopic subtypes reflect tumorigen- 
esis from distinct precursor cells, or represent distinct diseases that 
affect the same precursor cells via different environmental or 
epidemiological factors remains unknown. Nevertheless, our 
analyses reveal that they are sufliciently distinct to confirm the 
molecular heterogeneity of colorectal tumorigenesis and to merit 
consideration in cfinical management. 
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Table 1. Clinlcopathological features of colorectal neoplasms. 







All tumors 


HGD/T1 cancer 






/V=158 


/V=77 


Gender 


Male 


96 (63%) 


45 (58%) 




Female 


57 (37%) 


32 (42%) 


Age 


(yrs) 


68.1 


68.2 




(range) 


(37-89) 


(43-89) 


Tumor location 


Proximal 


75 (47%) 


35 (45%) 




Distal 


83 (53%) 


42 (55%) 


Size 


(mm) 


18.9 


22.4 




(range) 


(3-73) 


(5-50) 


Histology 


Adenoma 


118 


NA 




LGD 


81 (51%) 






HGD 


37 (23%) 






Tl cancer* 


40 (25%) 




Macroscopic type 


PN 


56 (35%) 


26 (34%) 




LST-G 


25 (16%) 


9 (12%) 




LST-NG 


48 (30%) 


21 (27%) 




S-FN 


10 (6%) 


2 (3%) 




DN 


19 (12%) 


1 9 (25%) 



*, all cases were submucosal cancers. Proximal, cecum, ascending and transverse colon; distal, descending and sigmoid colon, and rectum; LGD, low grade dysplasia; 
HGD, high grade dysplasia; PN, polypoid neoplasm; LST-G, granular type laterally spreading tumor; LST-NG, non-granular type LST; S-FN, small flat-elevated neoplasm; 
DN, depressed neoplasm; NA, not applicable. 
doi:1 0.1 371 /journal.pone.Ol 03822.t001 



Previous studies reported flat and depressed (refer to non- 
polypoid) neoplasms are associated with high risk of malignant 
potential [8,32]. However, in another study, no significant 
differences in the malignant potential between flat and polypoid 
adenomas were observed [33]. This discrepancy between studies 
could be attributed to differences in the macroscopic morphology 
of NPNs. In the present study, we showed the heterogeneity of 
molecular features among NPNs. Although the sample size was 
small, S-FNs showed few molecular changes associated with 
malignant progression, as reported previously by Morita et al. 
[34]. The remaining subtypes of NPNs had significant molecular 
signatures. Although a molecular analysis of NPNs was performed 
in several studies, different results were reported regarding 
molecular changes [17,3.5-37]. Clinically, DNs rapidly progress 
to invasive cancers even when they are small [38]. LST-Gs with 
large nodules or LST-NGs with depression tend to show 
histological invasion into the submucosa [39]. However, S-FNs 
are often characterized by the presence of low-grade adenomas. 
Therefore, morjjhological differences between NPNs should be 
considered in the design of appropriate screening and therapeutic 
intervention programs. Whether S-FNs could progress to other 
subtypes of CRNs remains to be determined. 

Recent studies have proposed a model for colorectal tumori- 
genesis that is consistent with three distinct molecular pathways 
(serrated, alternative, and traditional pathways) based on CIMP 
and MSI status, and BRAF/KRAS mutations [3,26]. We showed 
that macroscopic appearance was significantiy associated with the 
molecular phenotypes. Leggett and Whitehall [3] described an 
alternate pathway that is likely characterized by CIMP, KRAS 
mutation, and MSS, and they showed that no precursor lesion has 
been associated with this subgroup. However, the molecular 
features of this subgroup correspond to those of LST-G cases. On 
the other hands, DNs were characterized by LINE-1 hypomethy- 



lation and TP53 mutation, but not KRAS/BRAF mutation. 
CRCs with these molecular features commonly have a poorer 
prognosis than other CRCs [40] . This may be associated with the 
aggressive behavior of DNs. This macroscopic classification is also 
clinically relevant, since the macroscopic appearance of CRNs 
could be used to predict their molecular features. This could lead 
to better understanding of the pathogenesis of CRCs and improve 
the management of premalignaiit lesions specific to each 
macroscopic subtype. 

In the present study, we evaluated four genetic and nine 
epigenetic alterations. However, which molecular markers can 
identify the specific phenotypes of CRCs remains controversial. 
We selected markers to represent the distinct molecular features of 
CRCs. Alterations of these markers have been shown to be critical 
events during colorectal carcinogenesis [3]. However, few 
significant molecular abnormalities were observed in S-FN and 
LST-NG lesions using these molecular markers. Important 
challenges associated with the analysis of CRNs are heterogeneity 
between tumors, and poor reproducibility. Importandy, some 
studies may have been limited by bias caused by the small sample 
size; therefore, a large sample size is critical to examine the 
molecular features with adequate statistical power. In addition, 
possible disparities in the samples and sample collection methods 
among different studies, such as differences in the ethnicity of 
populations, and presence or absence of pretreatment may 
contribute to differing results. Advanced rectal cancer is treated 
with preoperative radiation in clinical practice, which could also 
cause bias and artifacts. Analysis of our study population shows 
that the sample size may be small and that there are different 
numbers of each macroscopic subtype. The small sample size 
could affect the estimate of the prevalence of molecular alterations. 
In addition, the smaller subgroups such as S-FNs and DNs could 
affect the power to detect specific molecular features between 
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Table 2. Cllnlcopathological and molecular characteristics among different macroscopic subtypes of CRNs. 









PN (%) 


LST-G (%) 


LST-NG (%) 


DN (%) 


P value* 






(N = 56) 


(N = 2S) 


(N = 48) 


(N = 19) 




Gender 


Male 


35 (63) 


12 (48) 


33 (69) 


12 (63) 


0.3873 




Female 


21 (37) 


13 (52) 


15 (31) 


7 (37) 




Age 


mean, yrs 


67.2 


71.2 


68.4 


68.4 


0.5900 




(range, yrs) 


(40-88) 


(55-89) 


(53-82) 


(43-87) 




Location 


Proximal 


22 (39) 


13 (52) 


25 (52) 


9 (47) 


0.5557 




Distal 


34 (61) 


12 (48) 


23 (48) 


10 (53) 




Size 


mean, mm 


16.9 


29.0 


18.8 


18.5 


0.0004 




(range, mm) 


(4-40) 


(12-73) 


(10-39) 


(7-35) 




Histology 


LGD 


30 (54) 


16 (64) 


27 (56) 


0 


<0.0001 




HGD + Tl cancer 


26 (46) 


9(36) 


21 (44) 


19 (100) 




Frequency of gene alteration/phenotype 


KRAS 


Mut + 


20 (36) 


17 (68) 


3 (6) 


3 (16) 


<0.0001 




Mut - 


36 (64) 


8 (32) 


45 (94) 


1 6 (84) 




BRAF 


Mut + 


1 (2) 


0 


0 


2(11) 


0.0378 




Mut - 


55 (98) 


25 (100) 


48 (100) 


17 (89) 




TP53 


Mut + 


3 (5) 


0 


6 (13) 


6(32) 


0.0028 




Mut - 


53 (95) 


25 (100) 


42 (87) 


1 3 (68) 




PIK3CA 


Mut + 


0 


4(17) 


1 (2) 


0 


0.0012 




Mut - 


56 (100) 


19 (83) 


44 (98) 


17 (100) 




MSI-H 


presence 


0 


0 


3 (6) 


2 (11) 


0.0713 




absence 


56 (100) 


25 (100) 


45 (94) 


1 7 (89) 




CIMP 


presence 


3 (5) 


8 (32) 


3 (6) 


3 (16) 


0.0028 




absence 


53 (95) 


17 (68) 


45 (94) 


1 6 (84) 




DNA methylation density (%) 


MGMT 


Mean 


13.0 


8.5 


8.7 


8.2 


0.0824 




95% CI 


9.0-17.1 


3.8-13.3 


5.0-12.4 


0.9-15.5 




SFRP1 


Mean 


49.5 


59.7 


39.5 


39.7 


<0.0001 




95% CI 


45.2-53.8 


54.7-64.7 


35.4-43.7 


34.8-44.6 




RASSFl 


Mean 


7.4 


6.4 


5.0 


5.4 


0.2562 




95% CI 


5.1-9.7 


3.5-9.4 


2.9-7.1 


1 .5-9.3 




LINE-1 


Mean 


63.2 


60.5 


61.4 


58.3 


0.0002 




95% CI 


62.0-64.4 


59.0-62.1 


60.2-62.6 


55.7-60.9 





*, P values were calculated by Chi-square test or Kruskal-Wallis test. PN, polypoid neoplasm; LST-G, granular type laterally spreading tumor; LST-NG, non-granular type 
LST; DN, depressed neoplasm; proximal, cecum, ascending and transverse colon; distal, descending and sigmoid colon, and rectum; LGD, low grade dysplasia; HGD, high 
grade dysplasia; MSI-H, high frequency microsatellite instability; CIMP, CpG island methylator phenotype; Mut+, presence of mutation; Mut-, absence of mutation. 
doi:1 0.1 371 /journal.pone.Ol 03822.t002 



subgroups. Our fmdings should therefore be confirmed in 
additional sample sets. 

Another limitation is that SAs were not included in this analysis. 
Many investigators have reported the relevance of SAs in 
colorectal oncogenesis [3-5]. Although both SAs, especially 
SSA/Ps, and non-SAs (conventional adenomas) frequently show 
flat-elevated morphology, they show different molecular features. 
Thus, we focused on the relationship between the macroscopic 
types and molecular features of conventional CRNs (non-serrated 
histology). However, we attempted to evaluate the clinicopatho- 
logical and molecular features in SAs on the basis of our previously 
reported data (Table S5 in File SI) [41]. Although the sample size 
was small, the frequency of CIMP tended to be higher in LST type 
than other types of SAs. However, we found no significant 
differences in any molecular alterations such as KRAS/BRAF 



mutations, and MSI-H among macroscopic subtypes. Flat- 
elevated type SAs were more frequently located in the proximal 
colon and LST-NG type SAs often showed the histological finding 
of SSA/P. 

Recent studies suggested that CRCs with KRAS mutations 
could be associated with a unique DNA methylation profile. 
CIMP2 CRCs exhibit increased age-related DNA methylation. 
Shen et al. [26] reported that the CIMP2 subgroup was highly 
correlated with KRAS mutations, but not MSI. Consistent with 
this, Ogino et al. [42] reported that CIMP-low (CIMP-L) tumors, 
which exhibit DNA methylation at a reduced number of CIMP- 
related loci, are significantly associated with KRAS mutations. 
CIMP-L tumors are also more frequentiy observed in men than in 
women. Hinoue et al. [43] identified the CIMP-L subgroup by a 
genome-scale approach and observed a higher frequency of 
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Table 3. Molecular alterations in 


relation to clinicopathological findings (multivariate analysis). 






Multivariate analysis 


Risl( factor 


Odds ratio 


95% CI 


lvalue 


KRAS mutation 


PN or LST-G 


9.11 


3.46-24.0 


<0.001 


Size (mm) 


1.07 


1.03-1.12 


0.001 


TP53 mutation 


LST-NG or DN 


5.30 


1.41-19.99 


0.014 


PIK3CA mutation 


LST-G 


26.53 


2.81-250.11 


0.004 


CIMP 


Age (yrs) 


1.14 


1.06-1.23 


0.001 


Size (mm) 


1.12 


1.05-1.19 


0.001 


LINE-1 hypomethylation* 


LST-G or DN 


3.41 


1.54-7.58 


0.003 


Histology (T1 cancer) 


4.40 


1.93-10.04 


<0.001 





*, We used the median of methylation density of LST-G and DNs (59%) as a cut-off value for LINE-1 hypomethylation. CIMP, CpG island methylator phenotype; PN, 
polypoid neoplasm; LST-G, granular type laterally spreading tumor; LST-NG, non-granular type LST; DN, depressed neoplasm. 
doi:1 0.1 371 /journal.pone.Ol 03822.t003 



KRAS mutations compared with other subtypes. However, the 
frequency of KRAS mutation was lower in CIMP-L subgroup than 
in CIMP2 (-50% vs. 92%). Yagi et al. [44] identified an 
intermediate-methylation epigenotype (IME) that was associated 
with KRAS mutations. They also demonstrated that CRCs with 
IME and KRAS mutation were correlated with poor prognosis. 
Although it remains unclear whether these classifications represent 
the same subgroups of CRCs, these subtypes appear to share 
overlapping characteristics. 

Global DNA hypomethylation may play an important role in 
genomic instability and colorectal carcinogenesis [45]. Several 
studies have shown that assessing LINE-1 methylation by 
quantitative pyrosequencing is highly reproducible, and the levels 
are correlated with global DNA methylation levels [40,46] . Here, 
we revealed that DNs are characterized by LINE-1 hypomethyla- 
tion and TP53 mutations. Consistent with these findings, 
Mastuzaki, et al. [47] showed that the LOH at 5q, 8p, or 17p 
was correlated with a lower level of LINE-1 methylation in CRCs. 
Furthermore, Ogino et al. [40] reported a significant association 
between LOH at 1 8q and a low level of LINE- 1 methylation in 
non-MSI-H CRCs and found that LINE-1 hypomethylation was 
independentiy associated with shorter survival times of CRC 
patients. In addition, TP53 mutations appear to be an important 
prognostic factor in patients with CRCs [48]. These findings 
suggest that the relationship between LINE-1 hypomethylation 
and TP53 mutation of DNs may be associated with their 
aggressive behavior. 

Yamauchi et al. [49] showed that the frequencies of CIMP, 
MSI, and BRAF mutation increase gradually along the bowel 
from the rectum to the ascending colon, suggesting the colorectal 
continuum. This phenomenon was confirmed by BRAF mutation 
analysis in an additional study by using a large-scale sample size 
[50] . In the present study, we evaluated the molecular findings of 
158 CRNs based on tumor location (Table S6 in File SI). 
However, the data corresponding to CIMP, BRAF, and MSI-H 
along bowel subsites was not consistent with the results of the 
above two studies. This may be due to difierences in tumor 
histology (adenomas vs. cancers or the exclusion of SAs) or 



different sample size between previous and present studies. 
However, we found that the frequency of KRAS mutation tends 
to be higher in cecal tumors than any other (75% vs. 14-37%), 
which is consist with the results of Yamauchi el al. [49]. In 
addition, 5 out of 6 cecal tumors with KRAS mutation were LST- 
Gs. These findings might be associated with site-specific tumor- 
igenesis. 

We have proposed the involvement of multiple parallel 
pathways in colorectal tumorigenesis (Figure 1) [26,51]. Briefly, 
age-related methylation such as SFRPl alters the physiology of 
colon stem cells, and predisposes them to acquire additional 
alterations. Subsequently, predisposed cells follow different path- 
ways on the basis of CIMP and MSI status, and genetic alterations 
such as BRAF, KRAS, and TP53 mutations. We observed 
increased methylation of SFRPl in all subtypes. However, our 
data revealed that LST-NGs and S-FNs showed no significant 
genetic or epigenetic changes. There are two possible explanations 
for the lack of characteristic molecular changes in these tumors. 
One possibility is that given that most CRCs develop through the 
adenoma-carcinoma sequence, some tumors could acquire alter- 
ations in cancer-related genes in the later stages of colorectal 
carcinogenesis. Yamamoto et al. [52] suggested that DNA copy 
number aberrations in CRNs occur as a late event in colorectal 
tumorigenesis. In addition, epigenetic alterations such as the 
promoter methylation of cancer-related genes begin to occur in 
precursor lesions such as adenomas. However, most promoter 
methylation events are more likely to occur during the transition 
from adenoma to carcinoma [53]. Another possibility is that 
tumors are more likely to involve other epigenetic alterations such 
as histone modifications [51] and microRNA changes [54], or the 
DNA methylation of other critical genes that may be involved 
earlier in the process of colorectal carcinogenesis. Additional novel 
molecular changes should therefore be investigated in LST-NG 
and S-FN cases. 

Several studies have demonstrated the potential of molecular 
biomarkers for the diagnosis or treatment of CRCs. Aspirin could 
suppress cancer cell growth and induce apoptosis by blocking the 
PI3K pathway. Importantiy, regular use of aspirin after diagnosis 
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Predisposed cell 
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Common 

CIMP/MSI Genetic or epigenetic 
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► 

Polypoid (protruding) / Cancer 



KRAS 
mutation 



(S> 



► 

Nonpolypoid (LST-G) / Cancer 



+ / - 



KRAS 
mutation 
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mutation 



LINE-1 
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(S) 
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LINE-1 
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Not determined 
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Figure 1 . Macroscopic subtypes in colorectal tumorigenesis. Precursor lesions can progress to cancer through the acquisition of epigenetic or 
genetic changes. Tumors from each subtype exhibit different characteristics, including their underlying molecular and genetic defects. However, 
whether small flat-elevated neoplasms can progress to other subtypes of CRNs remains unknown. PN, Polypoid neoplasm; LST-G, granular type 
laterally spreading tumor; LST-NG, non-granular type laterally spreading tumor; S-FN, small flat-elevated neoplasm; DN, depressed neoplasm. 
doi:1 0.1 371/journal.pone.Ol 03822.g001 



was associated with longer survival time in CRC patients with 
PIK3CA mutation, irrespective of aspirin use before diagnosis 
[55]. We observed a very low frequency oi PIK3CA mutations 
(3%) in our tumor samples (adenomas and Tl cancers), consistent 
with previous studies [37,56]. However, we observed an increased 
frequency of PIK3CA mutation in LST-G, compared with other 
types of CRNs, and four of these five lesions histologically revealed 
HGD or Tl cancer. Kaji et al. reported that PIK3CA mutations 
were detected only in LSTs with a higher pathological grade 
(Cancer or HGD) [56]. In addition, PIK3CA mutations are 
detected in 10-20% CRCs [28,57,58]. These observations support 
the hypothesis that LST-Gs may be premalignant CRC lesions 
with PIK3CA mutations and are associated with the PIK3CA 
mutation being a late observation in colorectal carcinogenesis. The 
regular use of aspirin may therefore suppress the transition from 
adenoma to carcinoma. 

In summary, we showed that CRNs could be divided into 5 
macroscopically distinct subtypes that differ in their DNA 
methylation status and genetic alterations, suggesting that different 
mechanisms are involved in colorectal tumorigenesis. However, 



further studies are required to clarify the epidemiology and clinical 
progression of the 5 CRN subtypes, which may have implications 
for the selection of optimal screening programs or therapies. 

Supporting Information 

Figure SI Endoscopic appearance of colorectal neoplasms (all 
lesions were observed after spraying with indigo carmine dye), (a) 
Polypoid neoplasm (O-I). (b) Small flat-elevated neoplasm (0-IIa). 
(c) Depressed neoplasm (0-IIc). (d) Granular type laterally 
spreading tumor (LST). (e) Non-granular type LST. 
(DOC) 

File SI Contains the files: Table SI. Summary of the 
macroscopic classification. Table S2. Clinicopathological and 
molecular features of small flat-elevated neoplasms. S-FN, small 
flat-elevated neoplasm; proximal, cecum, ascending and transverse 
colon; distal, descending and sigmoid colon, and rectum; LCD, 
low grade dysplasia; HGD, high grade dysplasia; MSI-H, high 
frequency microsateUite instability; GIMP, CpG island methylator 
phenotype; Mut-H, presence of mutation; Mut-, absence of 
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mutation. Table S3. The details of the multivariate logistic 
regressions, a) KRAS mutation: Logistic regression analysis using 
the stepwise method was performed to evaluate the relationship 
between KRAS mutation and gender (male vs. female), age (yrs), 
tumor location (proximal vs. distal), tumor size (mm), macroscopic 
type (polypoid neoplasm and granular type laterally spreading 
tumor vs. other types), and histology (Tl cancer vs. adenoma). SD, 
standard deviation; DF, degree of freedom, b) TP53 mutation: 
Logistic regression analysis using the stepwise method was 
performed to evaluate the relationship between TP53 mutation 
and gender (male vs. female), age (yrs), tumor location (proximal 
vs. distal), tumor size (mm), macroscopic type (non-granular type 
laterally spreading tumor and depressed neoplasm vs. other types), 
and histology (Tl cancer vs. adenoma), c) PIK3CA mutation: 
Logistic regression analysis using the stepwise method was 
performed to evaluate the relationship between PIK3CA mutation 
and gender (male vs. female), age (yrs), tumor location (proximal 
vs. distal), tumor size (mm), macroscopic t)'pc (granular t}pe 
laterally spreading tumor vs. other types), and histology (Tl cancer 
vs. adenoma), d) CIMP: Logistic regression analysis using the 
stepwise method was performed to evaluate the relationship 
between CIMP and gender (male vs. female), age (yrs), tumor 
location (proximal vs. distal), tumor size (mm), macroscopic type 
(granular type laterally spreading tumor vs. other types), and 
histology (Tl cancer vs. adenoma), e) LINE-1 hypomethylation: 
Logistic regression analysis using the stepwise method was 
performed to evaluate the relationship between LINE-1 hypo- 
methylation and gender (male vs. female), age (yrs), tumor location 
(proximal vs. distal), tumor size (mm), macroscopic type (granular 
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